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1.  Introduction 


The  aim  of  the  research  is  to  investigate  use  of  new  types  of  elastic 
waves  as  probes  for  examining  the  material  properties  of  piezoelectric  and 
electrostrictive  types  of  ferroelectric  ceramics  and  single  crystals.  This 
is  for  the  purpose  of  shedding  light  on  the  effectiveness  and  general  charac¬ 
teristics  of  fabrication  techniques,  as  well  as  exploring  basic  physical 
mechanisms  playing  a  role  in  the  technology  of  this  class  of  materials.  In 
the  category  of  relevant  physical  mechanisms,  one  can  cite  poling  phenomena 
and  grain  boundary  effects.  This  study  is  regarded  generally  as  influencing 
improvements  in  and  quality  control  of  existing  materials,  as  well  as  the 
evaluation  of  new  materials.  A  possible  spin-off  could  be  the  invention  of 
switchable  domain  wall  elastic  waveguides  for  electronic  signal  processing. 

During  the  present  report  period  our  effort  has  been  directed  toward 
developing  new  methods  for  analyzing  wave  guidance  or  scattering  from  a 
domain  wall  of  arbitrary  profile  in  a  material  with  general  elastic  aniso¬ 
tropy  and  piezoelectricity.  The  wave  guidance  (or  ducting)  of  acoustic 
waves  in  gases  and  liquids  has  already  been  extensively  analyzed,  using 
principally  the  WKB  approximation.  This  analytical  approach,  which  is  also 
widely  used  in  the  analysis  of  quantum  mechanical  energy  levels  for  an 
arbitrary  potential  function,  requires  that  the  mechanical  properties  be 
essentially  constant  over  distances  comparable  to  the  characteristic  wave¬ 
length  of  the  medium.  For  a  domain  wall  structure,  this  is  not  an  accept¬ 
able  limitation,  even  in  the  purely  elastic  case.  Furthermore,  the  WKB 
approach  turns  out  to  be  extremely  unwieldy  when  the  piezoelectric  effect 


We  have  developed,  for  general  symmetric  domain  structure  —  the  normal 
case  —  a  simple  numerical  trial  integration  method  that  treats  structures 
of  arbitrary  sharpness  and  anisotropy.  This  technique,  which  has  been  found 
to  converge  rapidly  and  inexpensively,  has  been  applied  to  both  the  iso¬ 
tropic  elastic  problem  and  to  the  simplest  known  case  of  a  piezoelectric 
domain  wall.  In  the  first  example  the  numerical  results  agree  with  WKB  when 
the  latter  technique  is  applicable,  and  in  the  second  example  the  results 
agree  with  published  results  for  wave  propagation  along  a  zero-thickness 
piezoelectric  domain  wall  of  the  same  type.  The  same  general  numerical 
trial  integration  technique  has  also  been  found  to  be  applicable  to  the 
calculation  of  wave  scattering  from  an  arbitrary  domain  wall  of  symmetric 
profile,  although  no  examples  have  yet  been  treated. 
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2.  Domain  Wall  Waveguiding 


2.1  Formulation  of  the  Waveguiding  Problem  for  a  "Duct"  in  an  Isotropic 


Medium 


This  year  we  developed  a  theory  to  handle  some  specific  cases  of  domain 
wall  waveguiding.  In  all  cases,  we  restricted  our  inquiries  to  shear  wave 
propagation  in  the  x-y  plane,  with  polarization  along  the  z-axis  and  axes 
chosen  such  that  the  domain  wall  region  was  centered  in  the  y-z  plane 
(Fig.  1).  The  regions  were  assumed  to  be  infinite  in  all  directions,  and 
the  fields  were  taken  to  be  uniform  in  the  z-direction. 

The  first  case  we  treated  was  waveguiding  in  a  nonhomogeneous  isotropic 
elastic  medium,  where  stiffness  c  (x)  and  density  p(x)  are  functions 

XJ 

providing  a  "slow  velocity"  region  (or  duct)  capable  of  guiding  a  shear  wave 
along  the  plane  x  =  0  .  This  problem  was  intended  to  be  only  an  overture  to 
the  more  pertinent  problem  of  the  piezoelectric  domain  wall.  Our  intent  in 
the  isotropic  case  was  to  model  a  domain  wall  as  a  region  of  lower  stiff¬ 
ness,  and  hence,  assuming  a  constant  density,  as  a  region  where  waves  propagate 
with  a  slower  shear  velocity  exists.  By  solving  the  acoustic  field  equations, 
we  hoped  to  gain  some  familiarity  with  the  type  of  modes  which  could  exist  in 
these  relatively  simple  isotropic  "domain  walls". 


By  writing  the  velocity  vector  as: 


where  v  **  v(x,y) 


itot 


and  assuming  that  all  fields  have  a  e 


time-dependence,  we  can  compute  the 


DOMAIN  WALL 


FIG.  1 — Basic  geometry  for  the  study  of  shear-type 
domain  wall  waves  in  an  infinite  medium 
uniform  in  7  . 


relevant  components  of  the  acoustic  strain  field  S  (where  J  -  1,2 . 6) 

J 

from  the  equation 


S 


J 


ico 


7  v 
Jk  k 


k  =  1,2,3 


(1) 


For  the  case  assumed  in  Fig.  1  we  have  S,  =  S  =  S.  =  S,  =  0  and 

12  3  6 


1  3v 
ico  9y 


—  S 


1  3v 
ico  3x 


(2) 


where  we  have  used  the  fact  that  3v/3z  =  0  .  The  stress  T  (I  =  1,2 . 6) 

can  now  be  written  using  the  constitutive  relation 


T 


I 


c  S 
IJ  J 


(3) 


where  c^  is  the  6x6  stiffness  matrix  and  each  component  is  a  function 
of  x  .  Using  our  expressions  for  in  this  equation  yields  1^=12= 
T3  =  T6  =  0  and 


T 


4 


c44  ico 


3v 

3y 


1  3v 
c44  ico  3x 


(4) 


Finally,  we  use  the  equation  of  motion 

Vkl  TI  =  1<0pvk  ’  (5) 

which  yields  two  trivial  equations  (0  =  0)  and 


3T  3T 

— —  +  — —  -  icopv  (6) 

3y  3x 
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Substituting  Eq.  (4)  into  Eq.  (6)  yields 


32v(x,y)  9  /  3v(x,y)\  2 

C44^  - 2 -  +  —  \C44^X^  - )  =  ““  v^x»y>  (7) 

3y  9x  \  9x  / 

-iBy 

Nov,  assuming  that  the  wave  propagates  with  an  e  variation  along  the 

y-axis  in  Fig.  1, 

v(x,y)  =  V (x)  e"i3y  (8) 

our  equation  becomes  an  ordinary  differential  equation  of  the  Sturm-Liouville 
type, 

£  (C*4(X)  ^Sl)  +  [<4(x)  ~  B2cA4(x)]  V(x)  =  0  .  (9) 

Equation  (9)  can  be  transformed  into  a  simple  Schrodinger  equation  by 
making  the  substitution 


ij>(x)  =  ^44<x)  V (x) 


g(x) 


2  p(x)  1 
co  - +  — 


c44(x)' 

2  1 

C44 (x) 

-C44(x>- 

2 

-C44(x>- 

and  now  becomes 


d  ^  2 

— 2  +  (g(x)  -  3  =  0 

dx 


(10) 


(ID 


which  can  be  solved  numerically,  or  under  the  right  conditions,  analytically, 
using  the  WKB  approximation.  We  will  not  report  the  results  for  this  case  in 
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detail,  since  most  of  the  isotropic  results  are  subsumed  in  the  piezoelectric 

case  which  follows.  Suffice  it  to  say  that  Eq.  (11)  is  a  completely  general 

result  independent  of  the  form  of  c^(x)  or  p(x).  Waveguiding  phenomena 

can  be  studied  by  assuming  appropriate  functional  forms  for  c^(x)  and 

p(x)  ,  so  as  to  provide  a  "duct"  at  the  plane  x  =  0  , 

As  a  general  qualitative  rule,  we  can  see  from  Eq.  (11)  that  waves  are 

2 

propagating  in  regions  where  g(x)  -  3  >  0  and  evanescent  in  regions  where 

2 

g(x)  -  3  <  0  .  This  follows  from  Eq.  (11),  and  the  fact  that 
2 

(x)  +  K  tKx)  =  0  has  solutions  of  the  form  'f'(x)  =  A  sin  Kx  +  B  cos  Kx  , 

it  2  ICx 

and  (x)  -  K  l p(x)  =  0  has  solutions  of  the  form  ip(x)  =  ce  +  De 

In  our  work  on  wave  guidance  (or  "ducting"),  we  want  to  have  a  function  g(»c) 

with  a  local  maximum  at  the  origin  (corresponding  to  a  shear  velocity  with  a 

local  minimum  at  the  origin).  We  see  that  there  will  be  a  maximum  and  a  mini- 

2 

mum  allowed  value  of  3  if  we  are  to  have  a  trapping  of  energy  (or  wave- 

guiding)  in  the  vicinity  of  x  =  0  . 

For  waveguiding  to  occur,  the  solution  to  Eq.  (11)  must  become  exponen- 

2 

tial  as  | x |  -*■  00  .  This  means  that  g(x)  -  3  <0  for  large  |x|  and  that 

2 

g(x)  -  3  =0  at  some  set  of  points  x  =  +  d  (Fig.  2(a)).  These  are  the 

turning  points  between  trapping  and  evanescent  regions.  Assuming  that 

Eq.  (11)  has  turning  points,  we  can  see  graphically  that  solutions  with 
2 

larger  3  are  damped  at  a  more  rapid  rate  as  x  increases,  and  have  turning 

points  which  are  closer  to  the  origin.  Therefore,  these  solutions  are  trapped 

2 

more  effectively  than  modes  having  smaller  3  .  Since  ft)  is  fixed  for  all, 

it  is  clear  that  the  slower  waves,  with  larger  3  ~  and  therefore,  smaller 
V  =  to/3  ,  are  the  ones  which  are  the  most  effectively  trapped. 


(c) 


FIG.  2 — (a)  Wave  is  trapped  in  region  |x]  ^  d  and  evanescent  elsewhere. 

(b)  Wave  cannot  propagate,  g(x)  -  3  <0  for  all  x  . 

(c)  Wave  radiates  laterally  (i.e.,  no  trapping),  g(x)  -  3^  >  0 
for  all  x  . 
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2.2  Waveguiding  at  a  Piezoelectric  Domain  Wall 


After  this  semiquantitative  preface  to  waveguiding,  we  now  include 
piezoelectricity  to  the  theory.  This  inclusion  is  necessary,  since  ferro¬ 
electric  domain  walls  necessarily  involve  the  piezoelectric  effect.  Once 
again,  we  assume  that  the  domain  wall  lies  in  the  y-z  plane  (Fig.  1) 
and  that  the  acoustic  waves  are  shear  waves  polarized  in  the  z-direction. 
We  now  study,  as  the  simplest  example,  a  hexagonal  6  mm  piezoelectric 
ceramic,  such  as  PZT,  with  its  poling  direction  taken  to  be  the  z-axis. 

The  constitutive  relations  can  be  written: 


T1  ~  C11S1  +  C12S2  +  C13S3  e31E3 


C12S1  +  °11S2  +  C13S3  "  C31E3 


C13S1  +  C13S2  +  C33S3  G33h3 


(12) 


and 


C44S4  "  C15E2 


C44S5  '  e15El 


C66S6 


e!5S5  +  G11E1 


e!5S4  +  £11E2 


(13) 


C31S1  +  e31S2  +  e33S3  +  C33E3 


where  J)  is  electric  displacement  and  £  is  electric  field.  We  now  wish 
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to  express  strain,  and  electric  field,  £  ,  in  terms  of  particle  dis¬ 
placement  14  and  electric  potential  <f>  .  Assuming  that  the  fields  are 
uniform  in  z  ,  we  write 


and 


a 


0 

0 

u 


where  u  =  u(x,y) 


4>  =  <Kx,y) 


Hence  S,  =  S  =  S  =  S,  =  0  and 
1  J  z  o 


9u  .  „  3u 

5  ~  3x  ’  4  3y 


(14) 


E 


1 


8<J> 

3x 


E3  =  0  (15) 


Therefore,  we  have  =  =  0  ,  and  our  only  nontrivial  equa¬ 

tion  of  motion  is 

3T  3T  32u 

— “  +  — “  =  P  — 2  (16) 

3x  3y  3t 

Also,  the  condition  V  •  J)  =  0  must  be  satisfied  throughout  the  region,  so 
that 

3D  3D 

—  +  — -  =  0  (17) 

3x  3y 

Substituting  Eqs.  (14)  and  (15)  into  the  constitutive  relations  (12)  and  (13), 
and  substituting  these  results  into  the  equations  of  motion  (16)  and  (17), 
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we  get 


a2 
9  u 

P  — 2 
3t 


0  = 


_9_ 

9x 


3 

9x 


E  3u  3<()  3 

C44  9x  +  el5  9x  9y 


E  9u 

C44  9y  ei5  9y 


(18) 


3u  S  3i 

15  9x  11  9x 


]  +  ¥ 


9u  S  9d> 
'15  9y  "  eil  9y 


Equations  (18)  are  simplified  by  taking  out  the  y  and  t  dependence 

u(x,y,t)  =  U(x)  e1^ 

<|)(x,y,t)  =  0(x)  e1 

E  S 

We  also  recall  that  p,  c^,  e^,  are  all  arbitrary  functions  of  x  . 

For  simplicity  in  what  is  to  follow,  we  will  define: 


C*VX) 


e^(x)  E  e(x) 


c(x)  ,  ei5^x^  =  e(x)  » 

The  equations  then  become  two  coupled  ordinary  differential  equations 


d 

dx 

dU 

C  dx 

J 

2 

•1-  (pco 

-  cp2)U 

=  - 

d 

dU 

-  e82U 

jd_ 

d4> 

dx 

e  dx 

dx 

dx 

_d_ 

dx 


d$ 

dx 


+  ee2<j> 


(19) 


-  e824> 


In  our  work,  we  wish  to  solve  these  equations  numerically.  By  making  the 
following  substitutions: 


dU 

dx 


=  W 


11 


we  obtain  four  coupled  first  order  equations, 


U'  =  W 

*’  =  0 

w'  =  $2U - j  |  epoj2U  +  (ee'+  ec')  W  +  (ee  -  ee'  )0 

c  e  +  e  ( 

9  1  (  o 

0f  =  3  $ - j  jepco  U  +  (ee*  +  cg’)  0  +  (ec*  -  ce*)  W 

ce  +  e  ( 

(20) 

where  the  primes  denote  differentiation  with  respect  to  x  . 

2.3  Numerical  Method  of  Solution 

We  want  a  general  technique  for  solving  Eqs.  (20).  The  technique  must 

2 

consist  of  choosing  an  appropriate  3  (the  eigenvalue)  such  that  a  trapped 
mode  does  exist,  and  then  determining  the  functions  U  and  $  from  numeri¬ 
cal  integration  of  Eqs.  (20).  In  our  piezoelectric  domain  wall  problem,  the 
only  change  in  the  ceramic  from  one  side  to  the  other  of  the  domain  transi¬ 
tion  region  is  a  gradual  reversal  in  polarization  P  .  If  we  can  describe 
this  polarization  profile  by  an  antisymmetric  function  P(x)  ,  then  we  can 
conclude  from  the  transformation  properties  of  the  permittivity,  piezo¬ 
electric,  and  stiffness  matrices  that  C(x)  and  c(x)  are  even  functions 
and  e(x)  is  an  odd  function  of  x  .  We  also  insist  that  the  solution 
functions  U(x)  and  4’(x)  be  smooth  and  continuous  throughout  the  region 
(for  a  smooth  and  continuous  polarization  profile  P(x)).  Moreover,  from 


12 


the  symmetries  and  antisymmetries  of  our  problem,  we  conclude  that  U(x) 

and  4>(x)  should  reflect  these  symmetries  and  antisymmetries.  Therefore, 

we  look  for  solutions  that  are  either  even  or  odd  in  x  . 

If  U(x)  is  even  (so  that  W(x)  =  U*  (x)  is  odd)  we  see  from  Eqs.  (20) 

that  4>(x)  must  be  odd.  Similarly,  if  U(x)  is  odd,  we  see  from  the 

equations  that  $(x)  will  be  even.  This  information  will  play  an  impor- 

2 

tant  part  in  determining  3  .  Going  back  to  Eq.  (20)  we  note  that  solu¬ 

tions  for  the  U(x)  even  mode  will  have  the  properties: 


U'  (0)  =  W(0)  =  0 

4>(0)  =  0 

at  x  =  0  and  solutions  for  the  U(x)  odd  mode  will  have: 


(21) 


«>'  (0)  =  0(0)  =  0 


U(0)  =  0 


(22) 


at  x  =  0  .  Our  numerical  integration  technique  consists  of  picking  trial 
values  for  U(x),  W(x),  $(x),  0(x)  at  a  point  z  =  a  >  0  that  is  suffi¬ 
ciently  outside  the  domain  wall,  in  a  region  where  analytical  solutions  can 
be  derived.  The  meaning  of  the  term  "sufficiently  outside"  is  that  the 

starting  point  should  lie  in  the  flat  part  of  the  profile  curves  in  Fig.  2. 

2 

With  these  "initial  values",  we  pick  a  value  for  3  ,  and  numerically 

integrate  our  equations  in  to  x  =  0  .  Note  that  the  initial  values,  as 

2 

well  as  the  actual  equations  of  motion,  depend  on  the.  choice  of  3  •  We 

2 

continue  to  pick  3  ,  each  time  integrating  to  x  --  0  ,  until  either 


Eq.  (21)  or  Eq.  (22)  is  satisfied  (depending  on  which  mode  symmetry  we  are 
seeking) . 

As  noted  above,  the  analytic  solutions  at  large  x  occur  in  regions 
where  c(x),  e(x),  e(x)  are  essentially  constant,  so  that  c*  (x)  =  e’(x)  = 
e'(x)  =  0  •  In  this  region,  Eq.  (20)  becomes 


U*  (x) 

(x) 


U(x) 


2 

j  epo> 

B  *(x) - 2  »(*) 

ce  +  e 


(a) 

(23) 

(b) 


These  equations  have  solutions  (for  x  large)  of  the  form 


where 


and 


U (x)  =  Au1(x)  +  Bu2(x) 
4>(x)  =  A<J>1(x)  +  B(f>2(x) 


(24) 


ux(x) 

u2(x) 

^(x) 

♦2(x) 


for  x  large,  positive 


— Kx 

Equation  (23(a))  clearly  has  only  one  solution,  u^  =  e  in  Eq.  (24), 
and  the  corresponding  potential  solution  is  related  to  u^  by 
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2  2 
substituting  <}»”  =  K  <J> ^  in  Eq.  (23(b)).  Substitution  of  U”  =  8  U  into 

— Bx 

Eq.  (23(a))  shows  that  a  solution  with  exponential  dependence  e  has 
U  =  0  ,  and  Eq.  (23(b))  then  gives  the  solution  =  e  in  Eq.  (24). 
From  Eq.  (24)  we  can  construct 


where 


W(x)  =  Aw^x)  +  Bw2(x) 
0(x)  =  A01(x)  +  B©2 (x)  , 

w^(x)  =  u^ (x) 

w2  (x)  =  u'2  (x) 

B^x)  =  ^(x) 

02(x)  =  <f>2(x) 


(25) 


In  Eqs.  (24)  and  (25)  the  constants  A  and  B  are  unknown,  but  their 
ratio  will  be  determined  by  the  conditions  in  Eqs.  (21)  and  (22)  at  x  =  0  . 
Since  solutions  to  a  set  of  linear  equations  can  always  be  multiplied  by  an 
arbitrary  constant,  only  this  ratio  B/A  is  significant.  We  note  that  the 
solutions  u^(x)  and  <f>^(x)  are  independent  of  the  solutions  u2(x)  and 
<J>2(x)  in  Eq.  (24),  and  hence,  we  expect  a  general  solution  of  the  form 


U(x) 

=  AU1(x) 

+  BU9(x) 

W(x)  = 

AWl(x) 

+  BW2(x) 

*(x) 

=  A«>1(x) 

+  B<I>  (x) 
l 

0(x)  = 

A01(x) 

+  B02(x) 
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for  all  x  >  0  .  (The  solutions  for  x  <  0  are  trivially  determined  by 
noting  the  condition  that  U(x)  is  to  be  even  or  odd.)  We  now  choose  a 
point  x  **  +a  ,  far  away  from  zero  and  take 

U^(a)  «  u^(a) 

ya)  =  4>^  (a) 

W1(a)  "  w^a) 

G^a)  =  01(a) 

defined  in  F.qs.  (24)  and  (25),  as  our  initial  values.  The  next  step  is  to 
numerically  integrate  Eqs.  (20)  in  to  x  =  0  ,  and  read  out  the  values 
VO),  ♦  (0),  W  (0),  0.^(0)  from  the  computer.  We  then  repeat  the  procedure 
for  the  functions  U2(x),  ^(x),  W2(x),  ©2(x)  ,  using  u2(a),  <J>2(a),  w2(a), 
©2 (a)  as  the  initial  values,  and  reading  off  U2(0),  $>2(0),  W2(0),  ©2(0) 
from  the  numerical  integration  to  x  =  0  . 

We  now  use  these  results  in  conjunction  with  Eqs.  (21)  and  (22)  to 
find  the  constants  A  and  B  .  For  U(x)-even  modes,  Eqs.  (21)  require 


W(0)  =  AW^O)  +  BW2(0)  =  0 
C>(0)  =  A4'1(0)  +  B4>2(0)  =  0 


(27) 


For  nonzero  constants  A  and  B  to  exist  and  satisfy  these  equations,  the 
determinant  of  coefficients  must  equal  zero,  i.e., 

W1(0)«t2(0)  -  W2(0)y0)  =  0  (28) 
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Similarly,  for  U(x)-odd  modes,  the  necessary  condition  for  a  solution  is 


©1(0)U2(0)  -  02(O)U1(O)  =  0  (29) 

2 

The  left  side  of  Eqs.  (28)  and  (29)  will  be  functions  of  3  ,  since  the 

integration  which  yields  Uj(0),  4>^(0),  U2(0)  ,  etc.  depends  on  the  value 
2 

of  3  chosen.  Hence,  we  define  the  functions 


F(32)  =  W1(0)$2(0)  -  W2 (0)^(0) 


(30) 


G(32)  H  01(O)U2(O)  -  02(O)U1(O) 

2  2  2 

The  values  of  3  which  zero  F(3  )  and  G(3  )  will  correspond  to  the 

allowed  even  and  odd  modes  which  exist  as  guided  waves  in  the  domain  wall. 

2  2 

In  general,  F ( 3  )  and  G ( 3  )  will  have  many  zeroes.  The  largest 
2  2 

3  which  zeroes  F(3  )  will  mark  the  lowest  order  U(x)-even  mode  allowed. 
2  2 

The  largest  3  which  zeroes  G(3  )  will  mark  the  lowest  order  U(x)-odd 
2 

mode.  As  3  decreases,  the  zeroes  alternate  back  and  forth  between 

2  2  2 
F(3  )  and  G(3  )  until  no  more  trapping  is  possible.  The  smallest  3 

which  zeroes  one  of  these  functions  will  be  the  highest  (and  most  weakly 

trapped)  mode  allowed.  By  this  procedure,  we  obtain  the  wave  numbers 

3q , 3j » 32 • • •  for  the  various  modes  propagating  in  the  domain  wall.  Once 

we  know  these  "eigenvalues"  we  can  go  back  and  read  out  the  values  of  U(x) 

and  $>(x)  by  numerical  integration. 


2 .4  Examples 

To  model  domain  walls,  we  must  be  able  tcj  choose  reasonable  representa¬ 
tions  for  c(x),  e(x),  and  e(x)  ,  the  stiffness,  piezoelectric,  and 
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permittivity  profiles.  It  stands  to  reason  that  these  functions  should  all 
depend  on  the  poling  profile  P(x)  in  some  fundamental  way.  This  dependence 
should  become  more  apparent  in  our  future  work.  For  now,  we  assume  that 
c(x)  and  e(x)  can  be  represented  by  Gaussian  like  profiles 


c(x)  =  %a(1  “  6e~“X  ) 
e(x)  =  e^l  ~  5e  aX  ) 


where  we  have  chosen  6  =  0.1  as  representative.  These  functions  have  a 
full-width,  half-maximum  given  by 


(32) 


We  will  refer  to  1  as  the  "width"  of  the  domain  wall.  The  function  e(x) 

is  also  a  function  with  a  width  f.  ,  and  is  intended  to  be  directly  pro- 

2 

portional  to  experimental  measurements  of  the  poling  profile  P(x)  .  For 
simplicity,  we  have  assumed  that  p(x)  is  constant. 

For  our  numerical  analysis,  we  have  chosen  PZT-8.  By  defining  dimen¬ 
sionless  coordinates  X  such  that 


(33) 


we  see  that  all  problems  treated  will  have  a  domain  wall  "width"  of 
AX  =  1  .  The  frequency  0)  will  be  "absorbed"  into  the  parameter  X  : 


(34) 
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For  each  problem,  we  specify  the  ratio  A/A  .  If  this  ratio  is  large,  then 
we  should  expect  strong  trapping;  if  it  is  small,  we  should  expect  weak 
trapping.  Finally,  we  need  to  specify  the  value 


R 


E  S 
C44  C11 
2 

e15 


(35) 


for  PZT-8,  R  =  2.3  (this  is  a  dimensionless  value).  Now,  we  can  see  from 
Eq.  (34)  that  for  a  fixed  P  ,  higher  frequencies  cu  result  in  larger 
A/A  and  hence  better  trapping. 

If  we  work  at  20  MHz,  then  A  858  100  pm  in  PZT-8,  and  we  see  from  the 
figures  that  we  can  get  fairly  good  waveguiding  all  the  way  down  to  about 
A  =  50  pm  (Fig.  3).  If  we  widen  our  domain  wall  to  A  =  100  pm  (Fig.  4), 
we  get  stronger  trapping  in  the  U(x)-even  mode.  We  also  note  that  a 
weakly  trapped  U(x)-odd  mode  appears  at  this  width  (Fig.  5).  Of  course, 
as  we  increase  A  ,  the  pre-existing  modes  become  much  more  strongly  con¬ 
fined,  and  more  higher  order  modes  appear  (Fig.  6).  For  A  =  200  pm  ,  we 
have  a  second  even  mode  and  a  second  odd  mode  appearing  (totaling  four  modes 
in  all).  Of  course,  if  we  were  given  a  fixed  domain  wall  of  width  A  ,  we 
could  reproduce  these  effects  by  doubling  frequency  (instead  of  doubling  A  ) 
to  get  better  trapping  and  higher  modes. 

The  next  question  which  we  address  is  whether  or  not  a  mode  exists  in 

the  limit  as  A/A  goes  to  zero.  Going  back  to  Eqs.  (30),  we  found  that 

2  2 
G(8  )  has  no  zeroes  for  A/A  <  1/2  and  F(8  )  has  only  one  zero.  This 

means  that  only  the  lowest  order  U(x)-even  mode  exists  as  the  domain  wall 

width  approaches  zero.  We  found  that  this  mode  persisted  all  the  way  down 

to  A/A  =  0.001  ,  at  which  point  the  numerical  program  accuracy  began  to 
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PHI(X)  U(X) 


PARTICLE  DISPLACEMENT  U(X)  VS.  X 


X 

ELECTRIC  P0TENTIAL  PHI(X)  VS.  X 


X 

FIG.  3 — Plot  of  l!(X)  vs  4>(X)  whore  the  abscissa  X  =  x/J.  is  expressed  in 
dimensionless  units,  and  where  the  ordinate  is  scaled  in  relative 
units  determined  by  the  normalization  U(3.0)  =  1,  4>(3.0)  =  1, 

and  H/X  =0.5  . 


particle:  displacement  urn  vs.  x 


X 


ELECTRIC  P0TENTIAL  PHHX)  VS.  X 
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FIG.  6 — Plots  of  U(X)  for  it/  A  ■  2  .  Two  u/en  modes  [(a),(c)]  and  two  odd 
modes  l(b),(d)]  exist. 


limit  the  reliability  of  the  result.  At  very  small  SL/X  ,  the  acoustic 
wave  sees  an  essentially  constant  c(x)  and  e(x)  ,  and  sees  a  piezo¬ 
electric  function  e(x)  which  looks  like  a  step  function  from  -e^,.  to 

H-e^,.  .  This  limiting  case  problem  has  an  analytical  solution  as  demon- 

3 

strated  in  the  paper  by  Maerf eld-Tournois .  The  solution  is  a  surface 

interface  wave  which  is  evanescent  in  x  and  propagating  in  y  .  The 
function  U(x)  is  even  and  <Kx)  is  odd,  which  corresponds  with  the 
assumptions  made  about  our  lowest  order  U(x)  even  mode.  Furthermore, 
it  can  be  shown  that  the  wave  number  3  predicted  by  Maerf eld-Tournois 
is 


where  X  and  R  are  defined  in  (34)  and  (35)  respectively.  This  result 

£ 

agrees  with  our  theory  as  y  -*■  0  ,  as  shown  in  Table  I. 


Table  I 


i/X 

Agreement  with 
Maerf eld-Tournois 

0.20 

4.1% 

0.10 

0.89% 

0.06 

0.12% 

0.04 

0.08% 

This  demonstrates  excellent  agreement  in  the  narrow  domain  limit.  Simi¬ 
larly,  U(x)  and  $>(x)  from  our  theory  show  excellent  agreement  with  the 
form  predicted  by  Maer  fc.l  d-Tournois  as  H/X  becomes  very  small.  We  con¬ 
clude  that  for  9JX  <0.1  ,  we  are  safe  in  dispensing  with  our  numerical 
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calculation,  and  using  the  simple  analytic  forms  given  by  Maerfeld  and 
Tournois: 


U(x)  =  Ae  ^ 


(37) 

« e  r  ~Kx  -3xi 

4>(x)  =  A  -  e  -  e 

for  x  >  0  ,  where 

2ir  1  3 

K  -  •  ■  i/2  _ 

\  [(R+l)(R+2))'  R+l 

and  3  is  given  in  Eq.  (36). 


2.5  Wave  Scattering  Problems 

The  procedure  developed  to  compute  guided  wave  modes  in  domain  walls 

can  be  also  applied  to  the  problem  of  reflection  and  transmission  of  waves 

propagating  through  the  domain  wall  (as  opposed  to  along  the  domain  wall). 

Of  course,  the  waves  propagating  outside  the  domain  wall  will  not  be 

evanescent,  so  that  new  initial  values  (i.e.  values  at  x  =  +a)  must  be 

2 

chosen.  Moreover  3  will  be  determined  by  the  equation 

3  =  k  sin  0 

o 


where  k  is  the  wave  number  far  outside  of  the  domain  wall  region,  and 
0q  is  the  angle  which  the  vector  Jc  makes  with  the  x-axis.  Therefore, 

3  will  be  smaller  than  before,  and  clearly  fixed  for  each  problem.  To 
use  our  technique,  we  must  restore  U(x)  symmetry  to  the  problem  so  that 
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we  can  use  Eq.  (21),  or  U(x)  antisymmetry  to  the  problem  so  that  we  can 
use  Eq.  (22).  This  is  done  for  symmetric  U(x)  by  launching  a  wave  of 
amplitude  A  from  x  =  -00  and  amplitude  A  from  x  =  +°°  . 

If  R  and  T  are  the  reflection  and  transmission  coefficient  respec¬ 
tively,  then,  in  the  region  x  >  0  ,  we  will  have  an  incoming  wave  of  ampli¬ 
tude  A  ,  and  an  outgoing  wave  of  amplitude  (R  +  T)A  .  The  sum  (R  +  T)  is 

then  chosen  such  that  Eq.  (21)  is  satisfied  after  integration  of  Eq.  (20) 
into  x  =  0  .  For  antisymmetric  U(x)  ,  we  launch  a  wave  of  amplitude  -A 

from  x  =  -00  and  amplitude  +A  from  x  =  +"  .  The  waves  in  x  >  0  will 

now  be  an  incoming  wave  of  amplitude  A  ,  and  an  outgoing  wave  of  ampli¬ 
tude  (R-T)A  .  The  difference  (R  -  T)  is  chosen  such  that  Eq .  (22)  is 
satisfied.  Hence,  from  (R+T)  and  {R-T)  we  can  uniquely  determine  R 
and  T  .  By  measuring  R  and  T  experimentally  at  different  angles,  we 
hope  to  be  able  to  measure  characteristics  of  the  domain  wall  itself.  This 
knowledge  could  enable  us  to  better  model  the  domain  wall,  and  could  be 
applied  directly  to  our  guided  wave  theory. 
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3.  Conclusion 


We  have  a  numerical  procedure  for  analysis  of  wave  guidance  or  scat¬ 
tering  at  domain  walls  of  practical  interest,  and  are  beginning  to  design 
some  experiments.  To  this  end  we  have  made  contact  with  the  ceramics  and 
materials  group  of  Professor  Eric  Cross  at  Pennsylvania  State  University 
and  have  agreed  to  work  co-operatively  on  the  basis  of  their  expertise  in 
materials  and  ours  in  elasticity.  Several  areas  we  wish  to  explore  together 
are: 

1)  Wave  guidance  experiments  in  the  new  strong  electrostrictive 
ceramics,  where  an  induced  piezoelectric  domain  wall  can  be  created  by 
applying  a  dc  electric  field  profile. 

2)  Wave  guidance  at  boundaries  in  layered  ceramics. 

3)  Possible  studies  of  macroscopic  models  of  grain  boundaries. 

We  hope  to  begin  with  (1)  immediately. 

4.  Presentations 

"Interaction  of  Acoustic  Waves  with  Ferroelectric  and  Ferroelastic 
Domain  Walls,"  American  Ceramic  Society,  Electronics  Fall  Meeting, 
Williamsburg,  VA,  September  16  -  19  —  to  be  prepared  for  publication. 

"Acoustic  Waves  and  Domain  Walls,"  Pennsylvania  State  University, 
University  Fark,  PA,  November  8,  1979. 

5.  Visits 

November  8  -  Pennsylvania  State  University,  University  Park,  PA. 
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